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ABSTRACT 


This thesis considers mathematical techniques for com- 
puting the optimal allocation of weapons from m different 
systems against n undefended targets. A standard nonlinear 
programming problem is considered. A discussion is given 
on John Danskin's Algorithm for the determination of the 
optimal values of the lagrange multipliers for this problem. 
Using a transformation of variables, the nonlinear problem 
is reformulated as a separable problem and solved by sepa- 
rable programming. A new method, the hybrid algorithm, for 
the determination of the optimal lagrange multipliers is 


deve Loped . 
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i ee NERODUCTION 


As a result of today's production costs, armament trea- 
[esc sand stockpiling facilities, both military tacticians 
and civilian defense planners are constantly confronted 
with the decision of how to optimally allocate offensive 
weapons such as ICBMs, SLBMs, bombs and artillery rounds to 
various military/industrial targets. Models with varying 
degrees of complexity and measures of effectiveness have 
been devised to assist in ‘Shi Geers t One 

This thesis concerns itself with solving these offensive 
weapon systems assignment problems so as to optimize measure- 
able returns such as damage or monetary savings. A non- 
immer model for an offensive allocation to a group of 
undefended targets is presented. An algorithm devised by 
John Danskin [Ref. 1] for obtaining the optimal lagrange 
fieipliers and hence the constrained optima to this probiem 
is discussed and illustrated with an example. The model is 
then transformed to a separable problem and an approximate 
solution is obtained utilizing separable programming. 

Using separable programming and the Kuhn-Tucker conditions, 
the thesis then proceeds to develop and apply a new method 
called the hybrid algorithm for calculating the optimal 


lagrange multipliers and then the optimal allocations. 
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ee oe MATHEMATICAL MODELS FOR THE ALLOCATION 
OF WEAPONS TO TARGETS 

The models considered in this paper are those that have 
a nonlinear objective function and linear constraints. The 
ilbmeec Specific models presented are those dealing with aerial 
bombing, "black and white" targets, and nuclear weapons. 
These three models have two similarities in common. First, 
the nonlinear function is directly dependent on the military/ 
industrial target value, the number of weapons allocated, 
and the weapon's Bere Givienless Oe ce Seveniesis against a 
meepet. secondly, the linear constraint restricts the 
number of weapons allocated so as not to exceed the total 
number stockpiled or available. | 

Similar weapon allocation models, considering some form 
of cost, have been developed. Examples of these are the 
allocation of weapons so as to inflict maximum damage with 
minimum delivery cost, or optimization of production costs 
subject to a budgetary constraint. This class of models 
will not be considered in this thesis, but any method 
presented may be modified to handle such cases. 

The basic assumptions relevant to all the following 
weapon allocation models are: 

1. The attack time is of short duration so as to 

preclude an effectiveness evaluation of preceding 


rounds. 





wee larcet. location 1s °fixed. 
Peeviilcaohe Kall by a single round 1s prohibited. 
4. The effectiveness of each individual weapon is 


Pic Mendel t. 


A. AERIAL BOMBING MODEL 

This model was developed by Koopman in Ref. 1. It is a 
model that might be used to allocate weapons, of differing 
magnitude, delivered by aerial bombardment. The total dam- 
mee intlicted by the bombardment is related to the number 
Seecdirect hits by a lethality function. This function re- 
@uees the target value by some fractional amount of its 
imevaous value. A lethality function representing y direct 


hits might be written as 


WGQ ie V jkr Wuciee O 1 Sil (1) 
Or in terms of damage as 
D(y) = V,(1-k*). 2 


If more than one type of weapon is used, the damage function 
takes on the form 


DO) = Vi (-fyk” 4). (3) 


Koopman assumes that each individual weapon (bomb) acts 
independently and that each bomb of type i has a probabil- 
ey Yas of hitting target j. Thus, the probability that y 


hits occur out of the x; bDOmbS dropped ony target } 1s 


“Vi 


P(y) = 7G 15 py%i-py) (4) 


and the expected damage takes the form 





Ae) MY NEGBA) |p. 


An application of Newton's binomial theorem reduces the 


expected damage formula to 


Jy, ]=¥ 


It 
— 
33 


Alen bAl| “ (6) 


Thus, the basic model to be used to allocate the aerial 
teapons 15 one that maximizes the expected damage subject 


to the number of bombs available and is stated as 


n i 
maximize 2, V Ls TT plik K 
pil isl 4 


subject to: 24, =X, where i=1,2,...,m (7) 
i 
—_ 
¥,= O and 
Ue 1S 5a sinee nec a. 
J 


In order to transform this model into a more mathematically 
tractable form the following approximation is utilized. 
Lemma If N is a large integer and  pC-¥)| is small, 


then 


Hl - /- - where =P y-2 





Proof 


eld] Fl ela 
pled} Lael 
me Ula 
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ee 1-| yk 


Hence, the basic aerial bombing model is now restated as 











maximize 5 v al exP|-LL. ; (8) 
Jed P Mh 
Serpgect to; 2.4. = x. Wine (C1 — lees... I 
z%. 20 
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B. "BLACK AND WHITE" TARGET MODEL 

This model was discussed by Danskin in Ref. 2 and is 
Seneilay toutnese publismed by denbrocder, Ellison, and 
Emerling [Ref. 3] and Mylander [Ref. 4]. A "black and 
White" target is one that is eitherdestroyed or not affected 
by the incoming weapon system, such as missile silos, or 
other undefended small point targets. Considering the case 
of m different weapon systems and n "black and white" 
targets, the kill probability of one weapon from each of 


these independent weapon systems is 





io (UEP, Yoo (la? Ne (9) 


Byeebecting, 


+f, = ma Shiga 
/-P=e le, 3 (10) 


(11) 


If there were xX; weapons allocated, expression (11) is re- 


wrrctLen as 
i= Beale: : —Bm¥m 
: CZ 


and the basic model for all the "black and white" targets 


would be stated as 
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C. NUCLEAR WARHEAD MODEL 

This model was first discussed and published by Lemus 
and David [Ref. 5]. It deals with the problem of allocat- 
ing nuclear weapons of varying yields to relatively defense- 
less targets. Associated with each weapon system 1s a 
launch reliability, Ge a nuclear yield, Ys, a circular 


error probable, CEP ., and a probability of penetrating 





th weapon system has been 


h 


defense j, Og Given the i 
successfully launched and penetrated and 3¢ defense,» tne 
target survival probability is 
ay) At" (ep) 
yrs 
if ee (14) 
where 
a and b are known constants 
ae is the hardness or measure of pressure which a target 
Gan resist. 
Hence the probability that a target will survive one attack- 
ing weapon 1s 
J -/-P-1-U R |r-4 
‘a 


ij : (15) 





By assuming a high launch reliability and relatively defense- 


ees target, it follows that 


a = 
at 


Cray 
anew '5) is restated as 
, P 
J=9 =A" 
a ee ole 


h 


Thus, the probability that the ie target will survive the 


attack is given by 


Spy. 3 
y, : | E; 2h he 
ey a 
=i J i (18) 
i 
where 47 [ A and Oi P 


10 





Hence, the basic nuclear weapon model, maximizing the value 
Of the targets destroyed, subject to the restriction that 
the total weapons allocated does not exceed the number 


available, is expressed as 


: 295 
maximize 2 Vio | (19) 
JES, J 
n 
Subject to: Er eee Glekee a IP en 
jel 
de 
ij 
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ee eee RMN Ad LON eOPeNECESSARY 
ANDVSULEICIENT CONDITIONS 

wie subsequent rok an calculations, deal only with 
models of the aerial bombing or "black and white" target 
form. If the reader desires to use the nuclear warhead mod- 
el, all formulas and programs must be changed to account 
for the new-constant X i. 

The necessary and sufficient conditions for the optimal 
solution to the weapon system allocation model, are obtained 
Ppaeaeadirect application of the Kuhn-Tucker theorem. The 
memcderm can find a complete discussion of the Kuhn-Tucker 


@emdations in Refs. 6 and 7. For the model, 


n 


maximize 2% V]reyp(2 yy.) 
imize a EP mi Pi i 


n 
Subject to: 2 Kw =. where 1=1,2,,...,m 
lee, 


XO 
1)j 
the Kuhn-Tucker theorem yields the following conditions 


* 
food 
x =O eee aeeeio | on acco) 
j'3 


BV, = 


Jj 


* * 
x > OSV, pe MiXj= A < HAV, (21) 
j j 
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* 
EnWs , > 1/H., ln M5Y; S510): 

J J Ne 
The following lemma shows, by contradiction, that the 
waomancge multiplier must be greater than zero and hence the 


eomstraint 1S active. 


emma the lasrange multiplier, {( Ye) = 0: 


Proof Assume } *=0. 
h POU. Ve 
then IN H; j 
thus x; 2 10)r 


buemiice LUN luteker conditions ,equlrTe that 


/ 


y= By Injt eo] = +6 >X XK. 


Hence, )\* must be greater than zero. The Kuhn- 
muGikey complementary Slackness conditions require 
that 

ye Xx, -X =O 


x 
Since * > Or then.2, xo, == “GED 
i J 








Extending the above to the model consisting of m weapon 
systems and n targets, the Kuhn-Tucker necessary and suffi- 


cient conditions become 


“f! oO ae V See Eily? ni" r (22) 
Z.. j OK, V jer Ey H x,[° A; E25) 


iS 





Momeany given target j}, there are three possible cases that 
wameoccur during the allocation process. Either the target 
will have no weapon system allocated to it, one weapon 

system allocated, or more than one weapon system allocated. 


x 
In the first case, for a given target j, ail ag and the 


Kuhn-Tucker conditions reduce to the form 


¥, OK, V =a Gl eONe ee eld S2lle? oo een (24) 


The second case requires that x, .=O for all weapon systems 


k 
kj 
fame... 1-1,itl,...,m and x52 0: Fourthissalloecatren: 


the Kuhn-Tucker conditions become 


HY ome LAy HN 








: (25) 
ees : |e AR 
PUP iS" PRM 
and on | See (26) 
Hy, | exp |p, | \ 
er -p ¥ |s oy 
PUP Ay Ay 
Hence . : * 27) 
Mid an 
ile. ay 
and ee ) In = (28) 


The final case is just an extension of the preceding one. 


For this case, the Kuhn-Tucker conditions become 


14 





ii Aafia a : % 
Herr amy x, for all i where a O (29) 


7 «| \* zt 
and Hy err aa for all k where ee (30) 
x 


: x 
Hence ———- = a constant for all x..>0. (315 
rij 7 





Which implies that 


x i 
EE Sg a icone metib al eS 0N x BOCs) 
Bij Ki a Kj 


. ’ 
and 29 My x, =) ea for all 1 where ar (35) 
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IV. NUMERICAL METHODS OF SOLUTION 


With the development of high speed computers many 
mimlyvercal and iterative methods have been devised for 
Wielding Exact and approximate solutions to the problem. 
Almost all methods treat the weapons as a continuous vari- 
able. Consequently an “eye ball" rounding approximation 
must be made for the final allocation. This thesis con- 
Siders the Danskin solution, which is an exact solution, 
the separable programming approximation, and the hybrid 
exact solution method. 

Other methods available, but not examined include 
Sequential Unconstrained Minimization Technique (SUMT) by 
Mylander [Ref. 4], geometeric programming by Passy [Ref.8], 
treatment as a transportation problem by Manne [Ref. 9], and 
Gener analytical algorithms by denBroeder, Ellison, and 


Emerling [Ref. 3], and Lemus and David [Ref. 5]. 


A. DANSKIN'S ALGORITHM 

one of the earliest algorithms for an exact solution 
was devised by John Danskin [Ref. 2] in the early 1950's. 
Danskin's algorithm uses the previously discussed Kuhn- 
Tucker necessary and sufficient conditions, Gibbs Lemma’, 
and the concept of marginal return, to obtain the optimal 
Ragrange multiplier. 


- Gibbs Lemma is discussed in Ref. 2. 
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The algorithm as applied to the one weapon system n 


target model consists of the following steps. 











Step 1 Comsi@erm tne Gduantlty ay and arrange the 
J 
t eumligs ie ane Vee Wels aN if 
arge IStine so a fa : aa I Oba 
a 
Seep 2 Deine stile saunG t1on 
A A j BY. 
s[AJ- Z\Al 3. In oe) 
pi HY? J x (34) 


where A denotes a trial value for the optimal 
% yas 
Voliicmeteticwl agranee mull Liplier nN arid cs aa 
J 
1s a trial allocation. Next find the largest 


index such that 
eee ee va) = (35) 
P5 J 
Denote this index as j=L. 
step 3 Onec Fos Knowmeand henceutmertarcets anethne 


x 
BOM LINd ean ecGl lS te ie » may be explicredy 


determined. From the Kuhn-Tucker conditions 


5\X}=X 








cee uV. 

X=2. In ale and trom (oC) 
ae r 

Sem 2 


A ie pV 
: Mall 2 Fy Bi In ie 2 


both of which contain the same list of regions. 





Subtracting (37) from (36) we obtain a formula 


for the optimal lagrange multiplier, 


ys 





-Lx-5 (1 ¥ 
ee =wV, i Ts 4 fl 
a (38) 


Kiowilno seme = index J, = may also be calculated 
by the formula 


L 
me HY )-x 


J= 


ce 


(39) 


Step 4 See mae has been found, the Kuhn-Tucker nec- 
essary and sufficient conditions can be 
itd ZedutO CDEAlieEne Opermal allocation. 

The algorithm as applied to more than one weapon system 
1s much more complex and requires a computer program to 
obtain the optimal target listing. A discussion of a 
iyerecal Computer program may be found in Ref. 10. Appendix 
Pemeome this thesis gives a heuristic approach to Danskin‘’s 


algorithm for one and two weapon system model. 


B. SEPARABLE PROGRAMMING 

Separable programming is used to obtain an approximate 
Folucron to nonlinear functions having a separable objec- 
tive function and constraint. A separable function 1s any 
general function that can be written in the form f (X)X,5-- 
= F(x), where £.(x,) HSoed LUC tom OF Cine Siig Le 


variable xX: A separable problem takes the form 


18 





n 
maximize 2. f.(x.) (40) 
1=l 1° 1 


n 
Bub ject to: ee ne) <b. OT alee ee 


ee 0 for all i. 


1 
The separable problem is then reduced to a linear program- 
fame problem by approximating each separable function by a 
peecewise linear function. References 6 and 7 should be 
consulted for a complete development of the piecewise lin- 
ear approximation. 

The separable programming formulation used in this 


thesis is commonly known as the lambda (A) method and 


takes the general form 


ijl ee 
maximize > SY \ £ (x) (41) 
j=l k=O kj kj jj 
: n rs < ‘ 
Subject to: = Jr g (x) b Ops l= 1..2 eel 
j=l k=0 kj ijk 3 1 
ee 
SJ » Bas | je Op Sales cad! pee amem mare (| 
k=0 kj 
ee 
k 


The lambda method Pree (minimizes) the actual piece- 
wise linear function values vice the slope formulation 
method that maximizes (minimizes) slopes of the approxi- 
Mating function. The slope method was tried for the weapon 
System allocation problem, but was abandoned because of the 
time and tedious effort required to update slopes for 
SUGCeSsive runs. 


eo 





The weapon system allocation problem (13) can be formu- 
lated as a aaa programming problem by introducing the 


new variable t.= 3 ml i -= problem (135) can then be restated 


i=] 32 
as¢ 
n -. 
maximize > V {/-e] (42) 
io 
Subject to: ue 2 Hy LO j=l. 2 eel 
a) 
4 Te fe Omi = Lat. 5 cee 
ay 
¥iy = O feorvalil- a Sandie 


minimize > VY ej (43) 
JJ 
subject to: t- “2 fOr falc. ee 
y* 
. 
2 4%-=X. Borer = 2 eee 
iy 
Ye =0Q for all ieee.) 


2 This approach was suggested by W.M. Raike, Assoc. 


Professor, Naval Postgraduate School. 
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By introducing a new variable Y, as 


os fal 








X 
“11 ned 


Pin Yon 


al 7Y 2n+1 


“aan Y(m+1)n, 
problem (43) is rewritten as 
A gi 
minimize > V, e~* 
k=| 


mM 


subject to: y 2, My Yak O 


S < 
Ai ni + Xx 


4.2 O 


for 


fOr 


The weapon system allocation problem is now in 


programming form 


men 
minimize > f 
k=! k 





4 


rea 


(44) 
k=l 2. Se 
Tu Il ae » Mm 


Ke dias le ele 


the separable 


(45) 


“4 





(m+)N 





Sepyeet Ol; 2 oO }-0 Oyen le see 
bi Se \ Jk 
(m+n 
DS | } =x for j=entl,n+2,...,n+m 
k=| ik Jk J 
ae 5 © 9 Res etd es EPR (20d 2 9 
where f,| | as Veo On Whe Wee 2 ose 
KK O (0) @a eons asa Amen wee Ona bN) 19 


VON) Gg) coat a eer 
Jk for k=j 


gu oel-4-Hjye for ten 


ainGleede= Meee. 2 oll 
O otherwise 
ane Or Gl tant 2 oes gl tm 
Y for k=n(j-n)+L 
= ; a k and £=1,2,...,n 
ik |Jk 
C) otherwise. 


Making the final transformation into the lambda form (41), 


the final problem to solve becomes 


minimize > > Auth Ly, (46) 


ZZ 





men "% 


Supeect £6: 
=| j= 


(m+n 


k=! i=0 


where 


tk fy 1 
O 


OMe Ag 2 5 ve sare 


Bi (Yk! 


O 
eid@er Ole | Tt... ean 
_) dk 
By (4 és 


i=0 \ Bip {y, |-O OTe Mier les 2 crt Th 


» A 5 $l = om ey tet ee ee tl i 


POM lal? ore. 5 (tn 


FOrewel 7 andes k 


iO ital. 4. 


for k=ntlant2. 34.5 (mem 
for k=} 

for k=ni+} 

Ae = SZ es 
otherwise 


POG k= (je 
Ame els 1. 2 eo al 


otherwise. 


emcee dll the separable, functions in the objective 


SoOuition. 


function of problem (46) are convex, the lambda separable 


programming method guarantees an approximate optimal 


Before discussing the method used to solve this problem, 


ZS 


it is necessary to understand the following definitions. 





Hetinition | Grid Size is number of increments into 


VittGmmugnenmamtCcnVal representing the range 


of ve is subdivided. 


teomimataom c Grid refinement 15 the process of increas- 


INGPtremOniaesS1 2c above  thit Of preeeding 


weeratrons. 


Definition 3 Nesting refinement is the process of re- 


GuGIng acne siailce sO Le tne sv aus ab ke Yr about 
Leamprcsent SOlmulon. | [ne erld size may 
increase, decrease, or remain constant. 
For this problem the range of en will be 
plus and minus one previous increment 


ADOMeGnENelMLeSele (Ss Ol Uta on., 


The optimal solution is obtained by the following lambda 


( ON ) method algorithn. 


oeep i Using the computer program contained herein, 


Beep 2 


generate the approximating function coefficients 
on punched cards suitable for the IBM MPS/360 
program, and a print out of the variable 

Gi gk KI FY foe Johto. ee nim. and k-n( pone 
for faces ANSPENeLtoOnsncOn proper Inpue 
dita OhtaiIned In the ecompucer program. 
Prac stiller oOu cpu sGards i homes cep lyin. bnew bn 
MPS/360 program. The A's associated with the 
855k") variable described in step 1 begin 
Vane Olumin nex (Grid Size +i) +r.  Sinee 


there is a one to one correspondence between 
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the A's and the Dr EMteoutatrom step 1. the 


solution for oF can be found by 


hOM- jag eae, Ml 
>< ’ > 
ik and i=1,...,m. 


= 9, 1V 
aT Bi 





Step 3 Repeat steps 1 and 2 using grid or nesting 
refinement, and continue until a minimum solu- 
trol sewOlbkarmeds srt should be noted that 
because or the plateau in the tail of the 
exponemt2al curve the nesting methods will net 
work if solutions lie in this area. 

In Appendix B, two numerical examples that were solved by 


this algorithm are given. 


C. HYBRID ALGORITHM 

Both Danskin's algorithm and the separable programming 
method have their advantages for the relatively small 
problems. However, for the larger problems these methods 
can become lengthy and cumbersome. To help alleviate this 
problem, the hybrid algorithm was developed. This method 
tii zes separable programming to obtain a trial target 
Mesteanad the Kuhn-Tucker conditions to calculate trial 
lagrange multipliers and solution. Refinements are then 
performed on this trial target list until the Kuhn-Tucker 
conditions yield optimal lagrange multipliers and solution. 

BeLore proceeding, it 1s necessary to define four 


merms peculiar to this algorithm. 


ZS 


Merinition 4 A link 175 a directed arc or branch connect- 


ing two weapon systems (nodes) via a shared 


ear 2 Cres, 


Meranition 5 A beta coefficient, 


/ fon. a= 1 


R o ce 7 G NEO Gaal # ean oti) 
i 4] au'nooes SP, «is link incident 
a oe to node i. 


Heranition 6 A Directed Tree 1S a connected graph which 


has no circuits or loops. Another defini- 
PROMS Gena da ¢raph @Setetree at and only 
if every pair of distinct nodes is connect- 


Ca Dy epmeGise lone path. 


Detamaition 7 A forest is a disconnected graph whose k 


CONMORCHES ake ELCes. 


The hybrid algorithm for solving the weapon system 


allocation problem, consists of the following steps. 


Seep | 


Secp Z 


Using the lambda separable programming method 
and a rather coarse prid, generate a trial 
Coeiialetancet la stinc wsa eM those targets 
where x55 20. Make a target by weapon system 
matrix with the element X denoting eo a and 
a blank for Se = QO. 

Draw all possible links between weapon systems, 
insuring that there are no more than one be- 
tween any two weapon systems and no weapon 
system has more than one link incident to it. 


This procedure will produce a forest of K 
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Seep 3 


Srep 4 


Beep 5 


N=ex p 


directed trees, where each weapon system repre- 
sents a node. The sum of the weapon systems 


in all the trees must equal the total available, 


K 
2m = m where Mm. om. (48) 


Decompose the forest and consider each tree 
separately. Designate a weapon system with no 
links incident to it as the base weapon system 
in each tree and relabel it number (1). 
Redesignate the other nodes (weapon systems) in 
the tree as 2,3,--+,M,. 
Calculate the beta coefficient (B.) for each 
weapon system in the tree, il.e. 1=1,2,...,m,. 
Calculate the optimal lagrange multiplier for 


the base weapon system from the following 













fe rm ba’; 
Me | / ; V We 
=| faliik k}— 
: aul Py) * iio Fi P 6 I=! is x 
Temes reek | 





Tar 






K 
i=! V¥,;208 i Pi kVSHARED /~i |Biy 
TETS NOT Lae TEM 
SHARED TGT 


(49) 
The above equation was derived by using induc- 
tive reasoning on various examples. It should 


be noted that the second summation in both the 
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: -_ 7 
oe 


=) wa 


se 





SEcp 6 


peep / 


pLcep 6 


Secp 9 


humerator and denominator is taken over all 
shared targets with the index 1 representing 
only one weapon system per shared target. 
Galculate ve for all the weapon systems in the 


tree from the equation 


* x 

Kad B (50) 
Revicatasteps 5» through 6 for all the trees an 
alc ea Oise Sta. 


Doanesthe miinm- Tucker conditions solve for ali 


* 
the ee bya the peollovine, Lommu lac. 





V 
* : 
V. pee In Ps where j 1S not 
J Bij Xr. shared (51) 
} 
and * | tal Whehe ) iS ssltaned 
fe alae and k is any system 
SHARING. JY eS oe - re 
Fer | sharing target j. 
: (52) 


Clicetmtorseceoitreal!) Tema mo wunne tucker Con 
CuelONsS ane —Satishied | lrasattsr led.) StOp. wan 
not, refine the separable program and repeat 


Steps sletnreuch 6. 


The hybrid method, like the separable programming 


method, has difficulty handling problems whose solution 
lies in the tail of the exponential curve. This is caused 
Pyaeene Computer's rounding errors and its inability to dis- 
Criminate between near zero values. The hybiwd. method: an 


its present configuration is incapable of solving problems 
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with lower bounds. However, by rederiving the Kuhn- 
moecker conditions and equation (49), this method can be 
modified to handle the additional constraints. Appendix C 


gives applications of the hybrid method to problems of vary- 


mr SiZeS. 


145) 


ce _ 
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V. EXTENSION OF MODEL AND SOLUTION 


The models discussed and used in this thesis are not 
inclusive nor necessarily representative of the present 
real world situation. Models are constantly changing in 
Mecder te reflect current technological and strategic 
@evamees. A model of current interest 1s one that considers 
an allocation of weapons against point and area defenses. 


Memea recent article addressing this problem, refer to 


ecmeort, f.A., Soland, R.M., Optimal Allocation of Missiles 
Sieest Area and Point Defenses ([Ref. 11]. Readers interested 


in building new models or modifying existing ones should 


fomsult Kooharian, A., Saber, N., and Young, H., A Force 


meecetiveness Model Wath Area Defense of Targets [Ref. 12], 


Perkins, F.M., Optimum Weapon Deployment For Nuclear Attack 


Meet 15), and Day, R.H., Allocating Weapons to Target 


Complexes By Means of Nonlinear Programming [Ref. 14], 
Eckler, A. and Burr A., Mathematical Models of Target Cover- 
age and tiissales Allocation [Ref. 15]. For a discussion of 


defensive models and methods of solution, the reader is 
referred to Dobbie, J.M., On The Allocation of Effort Among 
Meemrent Systems [Ref. 16], Brodheim, E., Herzer, I., Russ, 
mar, 48 General Dynamic Model For Air Defense [Ref. 17}, 


and Swenson, G.E., Anti-Ballistic Missile Allocations to 


Deeend largets With Time Varying Value Structures [Ref. 138]. 
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Whe reader interested in applying and extending the 
separable programming technique or the hybrid algorithm to 
other models, or models with additional constraints is re- 
Memnred to Ref. 19 and 20 for additional information on 


separable programming. 
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APPENDIX A 


The following two examples present an intuitive approach 
to Danskin's algorithm. Before proceeding with the examples, 
it is necessary to understand the concept of marginal return 
aioe inal utility). Marginal return is a measure of the 
Siemec in the objective function for a given change in the 
Baaependent variable. Symbolically this is expressed as, 

Maringal return = oS (A-1) 

meewe let X be or On A a continuous variable, 
then marginal return becomes, 

Marginal return = F! = ia (A-2) 
mimemerepresents the slope of the objective function at any 
given level x. It is assumed that any rational man will 
mimlecate nis independent variable so as to maximize (minimize) 


Micmieneinal return. With these ccncepts in mind, examples 


l and 2 are presented. 


WEAPON SYSTEM EFFECTIVENESS 


Example 1 










TARGET 


TARGET VALUE 









WEAPONS 
AVAILABLE 
Table 1 


Parameter Values 


SZ 





peep 


Seep Z 


[Witmeeteanlaro lide return ot the “Obj, ective 
MiMetlOmeWNinerespect co the independent 


variable oer ee. 


Marginal return = a = pV exp 


HH 





FOneenerduoverdatamtnlse ype las , 


MR (x, ) = -30exp(- .03x,) 
MR (x,) = -20exp(-. 2x.) 
MR (x2) = -10exp(-.02x,) 
MR (x) = -sexp(-.2x,) 
MR (xX, ) = ~2exp(-.2x,) 


Sollee earects wien tne Largest marginal retin 
aoe tile MOSt attractive and lucrative. the 
albocator will fire at those first until his 
memos are exhausted. ~ Whenealle the xis are 
zero, target one and two have the highest mar- 
(tire retians: 1se. “o0) did Z0sreSpeectively. 
Thus the weapons should be allocated to target 


imintil 1ts Marginal return equals that of 


Davee tied , 
-50exp(-.03x,) = -20exp(-.2(0)) = -20 
é 20 
exp(-.05x,) a ay 
xX, = WO Sal Se oes Or 
Since x —p15. 5109.50, more weapons can. be 
AelOea ted: 
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Now targets 1 and 2 are equally attractive 

Since they both have identical marginal returns, 
-20. Therefore, continue to fire at targets 1] 
and 2 until their marginal return equals the 
third highest, target 3, 


-30exp(~-05x,)=~20exp(~.2x,)=-10exp/-02(0}> -10 


xX, = eel si SO 102 
ay 7 eo lel 2) = 5.466 
DIMGe x era 40-0c0— 50, More weapons Can 


1 2 
be allocated. 


Now targets 1,2, and 3 are equally attractive 
Since they all have a marginal return of -10. 
inererores continue to-allocate to targets 1; 
2, and 3 until their marginal return equals 
the next largest, target 4, 


Nex t= .0SX-. i= -20exp(-.2x,)= -10exp(-.02x,)= 


1? 
-sexp|-.2(0) =e Ge 


Now XG io Be ead avaas| 


Bir 
Xz = SAO 
Since X, + X, + Xz = 101.314 >80, weapons can- 


ROUnDCM dlMbLOGated L£Omtaneetrs 1. 2 aana 5,5) Untied 
Eedean Mage tial We EUGTnMes Cauda 8tOn- 5. “lL here= 
fore, return to the point where their marginal 
return equaled -10 (called index L by Danskin), 


and find a place where their marginal returns 
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seep 5 


are equal and only 80 weapons have been allo- 
Giced ee liisewhtl Occur Lor a mMareinal return 
between -10 and -5S. 

et az be the Optimal marginal return between 

Wend o, samc sai tocate to targets 1,2 and 35 


itl et Neimoamdmomial Tetunrm equals =z. Thus, 








ZL; Oe =e ¥p ry OZ x) 





“Z=-seypl 08936 62 --zderp-2 




















/ / / = 
exp Osy |=exp\”.2%, = exp|—02%, |= 10 =K 


and hence x4 = eS ae 


1 -Ink 
X49 =e Vee 


X 2 = One 
Regardless of the constant K, where 2 ass ere ole 
used the same proportion of weapons x', x', and 
x. will be used to obtain any marginal return z. 
Therefore let K equal .9, and 
Rae af Se eV SEH? 


> 10997 J= .5268 


poles ROSA = 
x2 / 92 eens: 


Now find the proportions of weapons used. 
Seok 


Prepor 21 on Om x1 eee = oy 
: ; es ae 
Propervulon of x4 5-37 0566 
oe ae 
Puoportion-of x. > 37 a Oo 
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There are 80-40.086 = 39.914 weapons remaining 
to be allocated. Of the unallocated weapons 
37.7% belong to xj, 59.57% belong to x4, 


dice sor 05 Delong EO x Therefore the optimal 


t 
3 e 


Allocations are 


2K GOO eter oe 9 | 4p = = 541004 


il 


ra 
N 
It 


Sedbooee tt 2.0507 (59.9145 - 5.7245 


reed 5671592914) = 572.5 0ll 


TARGET WEAPON Ve volEv  EFPrPEGRIVENESS 
i 





TARGET 
VALUE 





5 ae? ae 10 
WEAPONS me) 
Rel be 
Table II 
Parameter Values 


Step 1 Calculate the marginal return of the objec- 


tive mUnet lol Wren nescpeet. te che Independent 


variable (x; 5). For the above data, this gives 
MR(x, 4) = -30exp -(..03x44 > .005x41)) 

MR (x, 5) = -20exp [- (2x1 = es 

MR(x,3) = -10exp |-"(.02x,, + 06x93) 


MR(x,4) = -Sexp [|-(.2x,,4 + 2x5 4) 


56 








2 





Be ade oe 





7 - a | — 
oe 
ak 
—_ Vie Tos 


i G 
+s 


ig 


_-* 


- 











etep Z 


MR(Xx,.) = -2exp F (.2x, = 2x, 6) 
MR(x,,) = “Sen COR, + .005x,, ) 
MR(x55) = -20exp - (2x4 1: Daa) 
MR(x, 3) = -30exp |- (.02x, , “7 06x, .)| 
MR(X54) = -~50exp vate =P oe 
MOS 3) = -2exp [-(.2x, , a 2x, .) 


Magecatteouvecanonesystem<lsto the S targets. 
iitemsahut21on ot example 1otves ene meededasolu- 
[MOMeNeamons are allocated CO targets 152, 
Hideominewl at herr Manecind! Gettin was —-0.500. 


The allocation was 


=-aSa, 
X44 ol. 67 6 
X19 = 53724 
Xz = 2225 Oe: 


Because of this allocation the marginal returns 


are adjusted as follows, 


MR(x, ,) = 6. 565exp(- .005x.,) 
MR (x, 5) = BS) eT 
MR(X, 3) = SOMES UGS) 
MR(x, 4) = ~Sexp(-2x,,) 

MR(X, ©) = -2exp(-.2x5-) 

MR(X, 4) = -1.06lexp(-.005x,,) 
MR (X45) = -6.565exp(-.2x,,) 
MR (X43) = -19.095exp(-. 06x, ,) 
MR(X, 4) = -50exp(-2x,,) 
MR(X5¢) = ~Zexp(-.2x,,)- 


Sf 


: ee ~ = 7 





Step 3 


meep 4 


PrOmecocmienToun tae target list allocating 
X, to tie targets with the largest marginal 
return, as previously discussed, until the 
problem is solved without a shared target or 
Mita lea Shared Garget WS Encountered... Hence 
allocate to targer 4 until 1ts marginal return 
equals the mext highest, target 3, 

- 50exp(-2x5,) = -19. 095exp |-.06(0)| = -19.095, 
Xoq = eh OT: 
There are 49.5203 weapons of type 2 unallocated, 
and target 3 has the next highest marginal 
BotMiimo Mee —taibvcteo adilkeadvanas weapons OL 
type 1 allocated, there is a possibility of a 
shared target. 
Target 3 will have weapons of type 2 allocated 
to it regardless of whether or not it receives 
weapons of type l. Tite rorone MEMmove bahee tao 
from xis target list and check to see if in 
fact it is shared. Adjusting the marginal 


returns to reflect the weapons of type 2 already 


allocated yields 


MR(x 1) = Urge Wee 

MR (x, 5) = -20exp(-.2x, 5) 

MR (Xx, 4) = -1.9095exp(-.2x, ,) 

MR(X) 0) = ~2exp(-.2x)_) 

MR(x,,) = -Sexp [-(.03x,, + .005x,,)| 
MR(x,,) = -20exp |-(.2x,, < os, 


Sy) 





MR(X>43) = -30exp |-(.02x,, ‘ 06X45) 
MR (X54) -19.095exp |-(.2x, ,) 
MR(X5 0) = -2exp |- (.2x, , + 2x, J]. 


Allocate weapon system 1 according to the newly 


If 


C2lculatcasnarcigal returms. “AS a result 
weapons will be allocated to targets 1 and 2 
Uilemimmrnei? Marginal returns equal -3.514, 

* = 71.492 


bw 


Xi * 8.694. 


Rdgustane the marginal returns to reflect the 


above allocation yields, 


MR (Xx, ,) = -3.514exp(-.005x,,) 
MR(x, 5) = ~3.514exp(-.2x,,) 
MR (x, 4) = ~1.9095exp(- 2x, ,) 
MR(x, 0) = -2exp(-.2x,,) 
MR(x,,) = =o Deter WD ocd 
MR(X,5) = ~3.514exp(-2x,,) 

MR (x, -) = es ES 
HSCs) = -19.095exp(-2x,,) 
ee = --Zexp(-.2x,,) 


Peep 5s Continue the allocation of weapon system 2 
bOmEene Larecet  t1St, VRecallethat target=4 nas 
-4797 weapons previously allocated. Allocate 
BOmEGhCCty os Until AfSs9margimaly return equals 
ao 0 95, 


-50exp(-.06x,,) = -19.095 


5 


S9 





x = 7.624 and 
IS 


x Pas ON 
24 os Z 
The marginal return for target 3 and weapon 


System 2 is now -6.55exp(-.0Zx. 2). Target 3 


aS 
Pose impossibly “shared sipee,-6.55 1S a) laroer 


Naan retUrn Ghan “5.5147 
Step 6 Continue allocating weapon system 2 until the 
Hempeinclerettirne or Lareet Sand 4 1s sequal to 
PneCmilcxGesmarOCSt. 2.) - 5. 5147 
-50exp(- 2x, 


J Se SYS G0) lope oer o ae 


4 23° 


thus 


X53 Son 120 


" 


Xo4 1.328 and 


X93 Se X54~ 42° 
The marginal return of weapon system 1 against 


target 3 1s now, 


13): 


See X Dp “| 02x, 


+ .06(35.726) =) il 76exp(-. 02x 


Simec the marginal return -17l7Gexp(-.02x._j< 


13 
oitexp(-.02X1-) , target 3 will not be shared 
in the optimal allocation. Continue allocating 
Meaponesystem 2 toethe target fist. larget 

Z, Which is the next largest marginal return has 
weapons of type 1 already allocated. Hence 
there is a possibility of a shared target. 


Seep /) Repeat steps 4 through 6 using target 2. Delete 


target 2 from weapon system 1's target list and 


40 


Seep 8 


allocate to remaining targets. Adjusting the 
marginal returns to reflect weapon system 2's 


Clb eieaadl LOcat on yieirus » 


MCS) = -s0exp|-(.03x, - 005%) 
MR (x, 3) = “1.17 6exp| = (023, + 06x.) 
MR(X, 4) = --S3514exp] (2x, + 24) 
MR(X) 5) = -2exp|-(. 2%) =: 26) 
MR(x,,) = -Sexp]-(.03% + 008x23)) 
MR(x,.) = -20exp |= (2), : 222) 

MR (x, 5) = “3.514exp] (02% es 06x)5) 
MR(x, 4) = “3.S1dexp| (2% 4 + 2x4) 
MR(X4,) - -2exp| (2%, + 2x, 4) ; 


As a result all 80 weapons can be allocated to 
@amcet 1 until 1ts marginal return equals 
-2.72exp(-.005x,,). 


Adjust the marginal returns to reflect this 


ClmTrente weapon system !allocatron, 


MR (x4) = -2.Texp |= (03%, + 005%7)) 


MR (x, 3) “1.176exp |-(.02%, + 06X93) 


MR =i. 5 ak ee + 2 
(x) 5 sex |= a 4) 


41 





Seep 9 


MR (x, <) = -2exp| (2%, + core) 


MR(x,,) = -.a83exp |-(.05x, ‘ 005x93) 
MR(x,.) = -20exp] = (.2X) + 22) 

MR (x, ) = -3.514exp] (02x, + .06x,.) 
MR (x, 4) = -3.514exp| (2%, ‘ 2x4) 
MR(X4¢) = -2exp| (2% : 2x, 4) : 


PhiI~Mocabe weapon sysuem 2 tO target 2 until its 
Niele Oniralleewe ul ane edlia l Saco olen AS 2 resi lt 
Xoo = SA? and X55 7 X43 a Xo4< X,- Thus, more 
weapons Gan be allocated to targets 2,3, and 4 
until their marginal return equals the next 


mancvest. -2,. Weapons ican ond, pe allocated 


Matil their Marginal returns. equal -2.8S6 


with 
X55 = 9,66 
oe = /59'-4 
Xo =e lie Ane 


Recalculate weapon system 1's marginal return 
for this current allocation. As a result 

MR (x, 5) = -2.8886exp(-.2x,,). 
pance Marginal return, Ob =-2-990 > -2./2, target 


2 Will be shared. Hence, the optimal solution 


for weapon system 1 is allocation to targets 1 
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and 2, and weapon system 2 is allocation to 
tareers 2, 5, and 4. 
* x 
sec? 10 The optimal marginal returns ONe and A>) can 
be calculated from equation (49), and the 


following solution can be obtained. 


WEAPON SYSTEM 


TARGETS 





Taipile sieri 
Weapon Allocation Solution 
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APPENDIX B 


The following two examples illustrate solutions by the 
separable programming method. Both solutions were obtained 


Meane the nesting refinement method. 


Eeample 1 


WEAPONS SYSTEM TARGET 
_ TARGETS EFFECTIVENESS: Mij VALUE 





Table IV 
ValwetoreParanebters 


VASES LS 202.259 : 20 2245.0 


15.0002 16205 : 16.0885 


G50 9 Ole si, : 6.435553 


Peooud bel. Zoe 





0 0 
Table V 
Weapon Assignment: X 


1j 
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Example 2 


WEAPON SYSTEM EFFECTIVENESS TARGET 
WARGETS it Z VALUE 





Table VI 
Value of Parameters 


79.14089 0 
.844646 Ser) O)ee1 2) 7 


39.5646 


LAZO So 
0 





Table VII 
Weapon Assignment After Four 
Runs Using Grid Size 20 
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APPENDIX C 


ime tollowing three examples illustrate solutions obtained 


by the hybrid algorithm. 


Example 1 








We SvCSM Msi Seales IME NSS TARGET 


VALUE 








TARGETS 


4 1 
5 5 1 S i 
mel a La 
AVAILABLE 1 i} 1 
Tapes yey 


Value of Parameters 


step 1 Use the separable program with a grid equal 
to .10, to obtain a target by weapon system 


Matrix. Draw all possible links. 


ice WEAPON SYSTEM 
TARGET 1 2 3 
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SEcpeaecaleculate all peta coefficients. 


fo =1 
Atl = 
Ale) 





* 
pine pu feca LCUlate Ne 


” LalS)+ ()ha(5) + In(é) + In) nl) if fats) ey sy ee 





=eX 
P a weet ae a pe a 
$ 5 5 a a 
* 
r | = «29544 


* 
ptcpeamecalculate all ROE 


ee Gee 

i Ne = .29544, 

SeChmommUoIny seme Kunn-lucker conditions calculace 
Bie .Opt iia leedtl LoOGat ion. and sinsume that. au 


ene Kuhn-fmcker Conditions ane Satistied. 


WEAPON SYSTEM ALLOCATION 
TARGET 1 2 5 





Table IX 
DOU Ono hat r 1x 
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Example Zz 


WEAPON OY SPEM EFFECTIVENESS TARGET 
TARGETS a Z 4 VALUE 





WEAPONS 
AVAILABLE | —20 


Table X 
Value of Parameters 





. 6 0 0 .48 0 
Bi 1 i not not 
_ shared |shared 


a pile sel 
Solution Matrix After 3 Separable 
Nesting Runs Using Grids of 
Sete Respectively. 
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Example 3 


WEAPON SY olEMPer PEC iiviEness TARGET 
TARGET 1 Z 3 4 > 6 VALUE 












10 we 156 ao 4 eo 5 900 
WEAPONS 
AVAILABLE; 10 20 15 5 10 > 


Table XII 
Value of Parameters 
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WEAPON SYSTEM ALLOCATION 





Tapie XL 
Solution Matrix After Separable 
RuneUstmo Gridesuvze 10. 


20 
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